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In this article, we study transportation network in Minnesota. We 
show that the system is characterized by Taylor's power law for fluctuation 
scaling with nontrivial values of the scaling exponent. We also show that 
the characteristic exponent does not unequivocally characterize a given 
road network, as it may differ within the same network if one takes into 
account location of observation points, season, period of day, or traffic 
intensity. The results are set against Taylor's fluctuation scaling in the 
Nagel-Schreckenberg cellular automaton model for traffic. It is shown that 
Taylor's law may serve, beside the fundamental diagram, as an indicator of 
different traffic phases (free flow, traffic jam etc.). 

PACS numbers: 89.75.-k, 89.75.Da, 05.40.-a 

1. Introduction 

In ecology, a striking observation is that variability in population abun- 
dance of a species and average population density are related in both space 
and time. The relationship is known as the Taylor's power law (or the law 
of the mean) and states that the mean (N) and the variance (A^) — (A)^ 
characterizing the number of population representatives are related by the 
power law 

(a2) - (A)2 = a{N)\ (1) 

with the characteristic exponent b describing effects of heterogeneity in spa- 
tial or temporal patterns of the frequency distribution. 

Since the publication of Taylor's famous paper in 1961 [Ij, the law of the 
mean has been substantiated in a vast body of ecological data that spans 
from protozooans to human populations [2] . This kind of fluctuation scaling 
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has been also noticed in a variety of natural phenomena including precipita- 
tion in a given area and river flows [3|. Taylor's law has also been observed 
in dynamics of different man-made systems driven by human activity [3]. 
Here, the list of documented examples includes behavior of the Internet and 
stock market dynamics. 

In this paper we perform a detailed analysis of fluctuation scaling in car 
traffic. We show that fluctuations in the number of cars passing through 
a given observation point fulfil the relationship given by Eq. ([T]). The clue 
observation reported in this paper is that the scaling exponent b does not 
unequivocally characterize a given road network. The parameter b may 
differ not only among various transportation networks but also within the 
same network if one takes into account: location of observation points (in- 
tersection, one-way street, highway, etc.), season, period of day, or simply 
traffic intensity. For this reason Taylor's law may serve, beside the funda- 
mental diagram, as an indicator of different traffic phases (free flow, traffic 
jam etc.). 

The outline of the paper is as follows: At the beginning (Section [2]) we 
analyze Taylor's fluctuation scaling in real data. In Section [3] we show that 
Taylor's law is also observed in the simplest cellular automaton model for 
traffic (the Nagel-Schreckenberg model), although interpretation of the ob- 
tained results deserves further studies. Section [J] gives concluding remarks. 

2. Fluctuation scaling in real traffic 

Hourly numbers of cars passing through observation points located on 
interstates, trunk highways, county state-aid highways, and municipal state- 
aid streets at various locations throughout Minnesota were retrieved from 
the Minnesota Department of Transportation [5]. The traffic intensity had 
been recorded by 72 automatic traffic recorders (ATRs) from 2002 to 2007. 
The datasets, from which mean and variance were calculated, include the 
number of cars observed by a single recorder and passing in only one direc- 
tion at a given hour in all weekdays of a month (consult Figure [T] for better 
understanding) . 

In Figure [51 one can see that although the whole data obey Taylor's 
fluctuation scaling with the characteristic exponent b = 1.43, in reality, 
the data are very heterogeneous and the parameter b characterizing traffic 
recorded by a single ATR may change dramatically from hour to hour. 
For example, the car flow (i.e. the number of cars counted within a given 
time window) as measured by ATR no. 222 possesses three distinct phases 
corresponding to non-overlapping periods of time by day, i.e. night hours: 
lam- 7am, working hours: 8am-6pm, evening: 7pm-12am. The phases have 
different values of the characteristic parameter b. It is meaningful that 
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ATR #222: 1am-2am; March 2004, direction North ATR #25: 2pm-3pm; February 2007, direction East 




Fig. 1. Two examples of time series considered in our analysis. The series represent 
number of cars counted by a single ATR during 20 working days in one month in 
the period 2002-2007. A) ATR no. 222, lam-2am, March 2004, direction: North. 
B) ATR no. 25, 2pm-3pm, February 2007, direction: East. 



during night hours the value of h = 1.12 is very close to unity 0, which 
suggest motion of independent (uncorrelated) cars. During working hours 
traffic is characterized by the largest value of 6 = 3.57. In this period, due to 
mutual interactions (correlations) between cars, fluctuations in traffic flow 
are signiflcant and congestion (traffic jam) can be observed. Our results 
clearly show that in the evening, nature of traffic recorded by the considered 
ATR changes again. The value of 6 = 2.36 attributing the new phase is lower 
than the corresponding value characterizing the earlier period. It implies 
that in the evening traffic becomes more homogeneous in comparison with 
working hours. 

Our results show that the characteristic exponent h in Taylor's power law 
for fluctuation may be used to characterize different phases of traffic flow. 
Larger values of the parameter characterize stronger collective phenomena 
manifesting themselves in congestion effects 0. We also show that Taylor's 
law characterizing various complex systems may comprise of several laws 
with different characteristic exponents. Doing such a decomposition one 
can better understand the considered systems. 



^ The scahng parameter 6=1 corresponds to Poissonian fluctuations (A^^) — {N)"^ ~ 
(N), cf. Eq. (m. 

^ It is worth to mention that a similar meaning of the parameter b has been also found 
in other complex systems. In ecology, larger values of b give evidence for stronger 
aggregation of individuals in different populations [T]. In genetics, they characterize 
stronger clustering of genes in chromosomes [6]. 



4 



Taylor'in'NaSch printed on September 10, 2009 



A 



B 



Transport 
in Minnesota 



Hours 



10^ r 
10^ r 



r 




■ lam -7 am 

A 8 am -6 pm 
• 7 pm - 12 am 




10-2 10-'' 10" 10^ 10^ 10^ 10"* 10^ 




mean 



mean 



Fig. 2. Taylor's fluctuation scaling in transportation network of Minnesota. A) 
Traffic intensity as measured by all automatic traffic recorders in Minnesota in 
2007. B) Daily fluctuations in traffic measured by a single recorder (ATR no. 222) 
in 2002-2007. In the figure, full points correspond to raw data (i.e. mean and 
variance calculated according to description given in the text), the open symbols 
express logarithmic binning of the data, and the solid lines represent their linear 
fits (in log- log scale). 

3. Fluctuation scaling in the Nagel-Schreckenberg model 



During the last 20 years, a number of cellular automaton models have 
been proposed in order to better understand complex traffic phenomena 
and reproduce the empirical data, such as the spontaneous formation of 
jams [ZllS]. A particularly simple, pioneering model for single-lane traffic 
has been proposed by Nagel and Schreckenberg in 1992 [9j. Further in the 
paper we show that the model, called the NaSch model, is characterized 
by Taylor's fluctuation scaling with non-trivial values of the characteristic 
exponent b. 

In the NaSch model a lane is represented by a one-dimensional array of 
L sites. Each site may either be occupied by one of N cars, or it may be 
empty {N < L). Each vehicle i = 1,2, . . . , N has an integer velocity Vi with 
values between and Vmax- State of the system at time t -|- 1 is obtained 
from its state at time t by applying the following rules to all cars at the 
same time: 



3.1. Description of the model 



(i) Acceleration: unless velocity Vi of a car is lower than Vmax it is ad- 
vanced by one 

Vi mm[vi + 1, Vmax]- (2) 
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(ii) Slowing down (due to other cars): if the distance di to the next car 
ahead is not larger than Vi the speed is reduced to avoid accidents 

Vi ^ mm[vi,di]. (3) 

(iii) Randomization: with probabihty p the velocity of a car (if greater 
than zero) is decreased by one 

Vi ^ Vi — 1 with probability p. (4) 

(iv) Car motion (update of positions): each car is advanced Vi sites. 

3.2. Simulation procedure 

Numerical simulations of the model have been performed on one-dimensional 
lattice of size L = 10^ with periodic boundary conditions. Each simulation 
has started at random initial conditions, i.e. random localization of N cars 
on a lane (further in the text we rather use the parameter p = N/ L than N) 
and random initial velocities. Then, we have updated state of the system in 
accordance with the rules described above. Collection of data has started 
after the first L time steps. 

3.3. Taylor's law in the model 

In our analysis we have concentrated on time series whose data points 
represent the number of cars passing through arbitrary lattice site in a given 
time period r (cf . Figure [1]) . Having such a time series one can calculate its 
mean and variance (both depending on r). Doing so for different window 
sizes r = 10, 20, 50, 100 . . . and plotting the obtained mean-variance points 
in a log-log graph we have obtained Taylor's power laws characterizing vari- 
ous stages of the traffic flow described by the set of model parameters Vmax-, 
p and p. In Figure [3] we have shown three scaling laws describing car motion 
in the NaSch model with Vmax = 10, p = 0.1 and different values of the car 
density p. 

In Figure UK and B one can see how the scaling parameter b character- 
izing fluctuations in the considered model depends on the car density p for 
two given values of the speed limit Vmax = 5 and 10 and the fixed probabil- 
ity of slowing down p = 0.1. Our simulations show that in the vicinity of 
the transition from the laminar traffic flow to the congested phase, it means 
in the vicinity of the maximum in the fundamental diagram (see Figure HP 
and D), values of the parameter b are significantly larger than outside the 
region. It means that the transition is turbulent. It is accompanied by large 
fluctuations, just like in continuous phase transitions. The sudden reduc- 
tion of the scaling exponent b exactly at the critical density (dashed lines 



6 



Taylor'in'NaSch printed on September 10, 2009 



1000 



100 

o 

« 

(5 
> 

10 

1 

1 10 100 1000 

mean 

Fig. 3. Fluctuation scaling in time series representing car motion in the NaSch 
model described by the parameters Vmax — 10, p = O.l and different values of the 
car density p. 

in the considered figure) is due to finite size effects, whicli cause tliat in the 
vicinity of the transition point traffic becomes locahy regular. 

The result are consistent with our previous observations obtained for real 
traffic, according to which below the transition point the growing density 
of cars causes higher values of b. These in turn show evidence of strongly 
inhomogeneous, turbulent traffic. 



4. Concluding remarks 

A short summary of the paper has been given in the Introduction (Sec- 
tion [1]), therefore we avoid it here. Nevertheless, we would like to draw the 
readers attention to issues related to the origin of the Taylor's law for fluctu- 
ation scaling. At present, universality of the law is widely recognized. The 
list of scientific disciplines in which Taylor's fluctuation scaling has been 
observed encompasses: genetics, epidemiology, medicine, physics, economy, 
computer and social sciences. It allows one to think, that there must ex- 
ist a common mechanism underlying this law. In our recent paper [10] we 
show that the conjecture is justified. In the paper, for the first time we give 
a universal, microscopic explanation of the law. Supported by real world 
observations ranging from insect and bird populations, through the human 
chromosome to traffic intensity in transportation network, we show that the 
law results from the density of state function characterizing the considered 
systems (the concept borrowed from statistical physics). 
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Fig. 4. Scaling exponent 6 as a function of the car density p in the NaSch model 
characterized by: A) v„iax = 5 and p = 0.1, B) v„iax ~ 10 and p = 0.1. Funda- 
mental diagrams describing efficiency of the car flow in the two considered cases: 
C) Vmax = 5 and p = 0.1, D) Vmax = 10 and p = 0.1. 
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